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The collisionless cold dark matter (CCDM) model predicts overly dense cores in dark matter halos 
and overly abundant subhalos. We show that the idea that CDM are decaying superheavy particles 
which produce ultra-high energy cosmic rays with energies beyond the Greisen-Zatsepin-Kuzmin 
cutoff may simultaneously solve the problem of subgalactic structure formation in CCDM model. 
In particular, the Kuzmin-Rubakov's decaying superheavy CDM model may give an explanation to 
the smallness of the cosmological constant and a new thought to the CDM experimental search. 



PACS numbers: 95.35.+d, 98.62.Gq, 98.70.Sa, 98.80.Cq 
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I. INTRODUCTION 

Recent cosmological observations such as dynamical 
mass, Type la supernovae, gravitational lensing, and 
cosmic microwave background anisotropies, concordantly 
predict a spatially flat universe containing a mixture of 
5% baryons, 25% cold dark matter (CDM), and 70% 
vacuum- like dark energy 0, 0] , termed as the standard 
ACDM model. The identities and the nature of dark 
matter and dark energy are among some of the biggest 
puzzles in contemporary physics. 

Although the nature of CDM is yet unknown, it is suc- 
cessfully treated in many aspects as weakly interacting 
particles. However, there exist serious discrepancies be- 
tween observations and numerical simulations of CDM 
halos in collisionless cold dark matter (CCDM) mod- 
els which predict too much power on small scales, 
manifested as cuspy CDM cores in dwarf galaxies 0, 
galaxies like the Milky Way 7], and central regions of 
galaxy clusters |8| as well as a large excess of CDM sub- 
halos or dwarf galaxies within the Local Group 0- 

To alleviate the discrepancies, among many other at- 
tempts, models of non-standard interacting CDM have 
been proposed. They include self-interactions 0, annihi- 
lations ^(j, and decaying cold dark matter (DCDM) 0, 
H^ . Although these models involve different interactions, 
almost all interactions result in an adiabatic expansion of 
the cuspy halo that lowers the core density and reduces 
the number of subhalos. However, both self-interacting 
and annihilating CDM models require embarrassing large 
interaction cross-sections that have made the models less 
appealing. Although DCDM models are viable, possible 
underlying particle physics has been ignored. 

Another big puzzle in astrophysics is the origin of the 
ultra-high energy cosmic rays (UHECR). One may ex- 
pect that UHECR should originate from some unknown 
astrophysical sources at extragalactic scales. Greisen, 
Zatsepin, and Kuzmin (GZK) observed that due 
to inverse Compton scatterings of the relic photons the 



UHECR energy spectrum produced at cosmological dis- 
tances should steepen abruptly at energy ~ 10 10 GeV. 
However, a number of cosmic ray events with energies be- 
yond the GZK cutoff have been observed by Fly's Eye fl4| 
and AGASA 0] . A simple solution to this impasse is to 
invoke new physics in which UHECR can be produced in 
a cosmologically local part of the Universe. Ideas such as 
long-lived metastable superheav y p articles that are de- 
caying at the present epoch [T3, OH uil 1 annihi- 
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rations of stable supermassive particles in halos [21| , and 
collapses of cosmic topological defects j22| have been pro- 
posed. In most of the models the superheavy objects can 
simultaneously be viable candidates for DM. 

In this paper, we try to address these issues at the 
same time within a single theoretical framework. We 
pursue the DCDM scenario, suggesting that the CDM 
is decaying weakly interacting superheavy particles with 
mass of the grand unification scale. In our scenario, not 
only the decay would produce much less concentrated 
cores in CDM halos, but also the decay products contain 
highly energetic quarks and leptons which lead to the 
production of ultra-high energy cosmic rays (UHECR) 
with energies beyond the Greisen-Zatsepin-Kuzmin cut- 
off. Moreover, the longevity of the superheavy particles 
may shed new light on the origin of the observed small 
value of the cosmological constant. 

The paper is organized as follows: In section II we 
illustrate our idea by using the Kuzmin-Rubakov model. 
After briefly reviewing this model, we show in section III 
how this model can be naturally fitted into the scenario 
of DCDM. We show how this model solves the cuspy halo 
problem, and find out the parameter space which allow 
us to solve the origin of UHECR as well. In section IV we 
discuss some phenomenological implications and suggest 
that some on-going experiments could test this scenario. 



II. KUZMIN-RUBAKOV MODEL 

Here we will concentrate on a specific scenario pro- 
posed by Kuzmin and Rubakov (KR) |19| and show how 
the KR scenario for producing UHECR is related to the 
subgalactic structure of the Universe. 

KR [13 considered an extended standard model with 
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a new SU(2)x gauge interaction and two left-handed 
SU(2)x fcrmionic doublets X and Y and four right- 
handed singlets. Here at least two doublets are intro- 
duced because the SU {2)x anomaly prevents the number 
of SU(2)x doublets from being odd. All new particles are 
singlets of the standard model, while some conventional 
quarks and leptons may carry non-trivial SU {2)x quan- 
tum numbers. The SU{2)x gauge symmetry is assumed 
to be broken at certain high energy scale, giving large 
masses rax.Y to all X and Y particles. Furthermore, 
X and Y are assumed to carry different global symme- 
tries, so there is no mixing between them. As such, both 
the lightest of X and the lightest of Y, which we call X 
and Y respectively, are perturbatively stable. However, 
SU(2)x instantons induce effective interactions violating 
global symmetries of X and Y. Assume mx > ffiy, then 
the instanton effects lead to the decay 

X — > Y + quarks + leptons (1) 

with a long lifetime roughly estimated as Tx ~ 
m x 1 e^l ax , where ax is the SU(2)x gauge coupling con- 
stant. With the choices m x > 10 13 GeV and a x J$ 0.1, 
tx ^ lOGyrs and X particles are decaying at the present 
epoch. There have been many discussions on the pro- 
duction of X particles in the early Universe. X particles 
may be produced thermally during reheating after in- 
flation with the produced energy density comparable to 
the critical energy density of the Universe [jjj (see also 
Refs. 0,Hf|). Also, it was realized in the same or differ- 
ent context that superheavy particles can be efficiently 
generated from vacuum quantum fluctuations during in- 
flation l2ol or couplings to the inflaton field during pre- 
heating |27j . 

The particles X and Y are good dark matter candi- 
dates. According to KR, there are two possible outcomes 
after X particles have decayed. If Y particles are pertur- 
batively stable, they are also stable against instanton- 
induced interactions in virtue of energy conservation and 
instanton selection rules. In addition, if mx ^ Tny, an 
approximately equal amount of Y particles is produced in 
the early Universe. Therefore, the decay products would 
contain stable supermassive Y particles that constitute 
a dominant fraction of the CDM with a small admix- 
ture of X particles as well as highly energetic quark jets 
and leptons that subsequently produce UHECR. Alterna- 
tively, the Higgs sector and its interactions with fcrmions 
may be organized in such a way that Y particles are 
in fact perturbatively unstable. As such, Y particles 
would instantly decay into relativistic particles and leave 
metastable X particles being the CDM. 

Intriguingly, it has been recently pointed out that if the 
longevity of the superheavy particles in the KR model is 
due to instanton-induced decays, the observed small but 
finite cosmological constant can be explained by instan- 
tons or vacuum tunnelling effects in a theory with degen- 
erate vacua . In such a theory, the vacuum energy 
density of the true ground state is smaller than that in 
one of the degenerate vacua where we live now by an ex- 



ponentially small amount if quantum tunnelling between 
the degenerate vacua is allowed [24j| . 



III. RESOLUTION OF THE CUSPY HALO 
PROBLEM AND UHECR 

We now turn to the cuspy halo problem and show how 
this problem can be solved within the context of the 
KR model. Numerical simulations of CCDM halos show 
cuspy halo density profiles well fit with the generalized 
Navarro-Frenk- White (NFW) form HHQ, 

+ (2) 

with the slope parameter a ~ 1 — 1.5 and the concen- 
tration parameter c = r2oo/r c ~ 20, where r c is the core 
radius, p c is the mean density of the Universe at the time 
the halo collapsed, and ^200 is the radius within which the 
mean density /?2oo is 200 times the present mean density 
of the Universe. However, observations indicate flat core 
density profiles with a < 0.5 and smaller concentrations 
with c ~ 6 — 8 0,0,0. Below we will simply study the ef- 
fect of DCDM to the original NFW profile with a = l@ 
in Eq. Defining x — r/r2oo, it gives the halo mass 
profile M(x) = M 2 qoF{x) that is the mass within x and 
the associated rotational velocity V{x) — V2oo[F (x) / x]* , 
where M200 = M(x = 1), V200 = V{x = 1), and 

Fix) = [ln(l+or)-ca;/(l+CT)]/[ln(l+c)-c/(l+c)]. (3) 

Suppose a CDM halo gas composed of X particles is 
formed at some high redshift with the NFW profile and 
a velocity dispersion vx = \J GMwo^x /2r2oo.x , where 
Maoo.x is the mass of X particles within the radius ?"2oo,x- 
The observed velocity dispersion typically ranges from 10 
to 1000 km/s for dwarf halo to cluster halo. In X's rest 
frame, the decay Q produces a Y with a recoiling veloc- 
ity 7 rc w rc = (5(1 - S/2)/(l - (5), where 7 rc = ~ v% c 
and S — {mx — my)/mx, and highly relativistic quarks 
and leptons of energy E q i = 7 rc w rc'™j s :(l — 5). The value 
of 5 depends on the detail dynamics of the high energy 
model. Here we will treat it as an input parameter. There 
are two possibilities. When 1 > 5 > vx, we find that Y 
would be relativistic and/or beyond the escape velocity 
of the halo. This together with the case of an unstable 
Y correspond to the scenario discussed in Ref. [ll|, to 
which readers may refer for details. In the following, we 
will discuss the case for S < Vx, i.e. nearly degener- 
ate masses, in which stable Y particles would be bound 
to the halo with an averaged velocity about \Jv' x + vf c 
(v rc ~ 8) just after the decay of X particles. In particu- 
lar, S ~ 1 — 2 x 10~ 4 corresponds to the case considered 
in Ref. \V2}. 

Let us assume that most X particles have decayed and 
that the halo of Y particles with the NFW profile has 
been formed by now. Using the virial theorem it can be 
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shown that the core radius has expanded to 

1-2J S 2 (l-S/2) 2 



Tc.y = r c ,x/y, 



y 



i 



X 



(1-5)= 



(4) 



We will follow the method in Ref. |llj to work out the 
consequences of this core expansion. The difference is 
that here the mass inside r2oo,x /y is only slightly changed 
to (1 — S)M2oo,x- As such, the final density within 
r20Q,x/y is y 3 {l - 5)p2oa- To obtain r 20 o,y, we solve for 
r = y r 200,Y hi Eq. P)l (a = 1) within which the initial 
density is y _3 (l — 5)~ 1 P2W- The resulting equation is 



,0.2, 



x 3 F^ 1 (x) = j/ 3 (l— S) and we find that r2oo,Y — y""f2oo,x 
for y < 1 and 6 << 1. Hence we obtain cy ~ y 12 cx- To 
circumvent the over-concentration problem, y should be 
about 0.4, implying that 5 ~ Q.77vx ■ Using cx = 20, y = 
0.4, and Eq. J3J), we obtain the mass profiles and rotation 
curves of the original X halo and the presently formed Y 
halo shown in Fig. 1. We find that M2oo,y — O.58M200.X 
and M Y (r = 0.1r 200 ,y) ^ 0.27M x {r = 0.1r 200 ,x), and 
that V 2 oo,Y ^ 0.83y 2 oo,X, Ki, B:t| y ~ 0MV max , x , and 
~ 2.5r max ^ x , where U ma:E is the maximum rota- 
tional velocity at radius r max . In Fig. 1, we have also 
reproduced the mass profile and the rotation curve for 
the case ^lj in which X decays into relativistic parti- 
cles. This requires solving for x = r/r2oo,x in the equa- 
tion x^F^ 1 (x) — y 4 , where the mass inside T2oo,x/y is 
yM2oo,x and y — 0.5 is the fraction of X particles that 
still remain by now. In this case, the softening of the 
central concentration is the same as in the Y halo, but 
the reduction in the halo mass profile and the flattening 
of the rotation curve are even more pronounced. Thus 
we have shown that one can put KR model which was 
originally proposed to explain the origin of UHECR into 
the DCDM model. 

Now let us examine the production of UHECR in the 
scenario proposed here and the applicability of the virial 
theorem for obtaining the Y halo profile in Eq. It was 
found that the level of the UHECR and the UHE neu- 
trino fluxes produced from X decays is proportional to 
a single parameter rx = £,xtol T x for a fixed mj, where 
£x is the present fraction of X particles in CDM and 
<o = 13.7Gyrs is the age of the Universe and there 
rx = 5 x l0~ n /5 was used to fit the observed UHECR 
flux spectrum. Note that a factor of S is added because 
the energy of the decay relativistic quarks and leptons is 
E q j ~ Smx, where S ~ 1 for the case in Ref.fl8| and here 
S ~ 0.77v x ~ 10~ 3 (where v x is about 300 km/s) and 
mx = 10 16 GeV, and also that the parameter rx will be 
larger if the energy dissipation of the decay particles is 
taken into account [2j| . Assume that the X halo is origi- 
nally formed at 0.1 — 1 Gyrs and that tx = 0.7Gyrs. Then 
the dynamical effect of X decays on the halo is at work 
from about 0.7 Gyrs to the present time. Since X and Y 
are non- relativistic and tx << to, most X particles have 
decayed into Y particles many Gyrs ago and the Y DM 
halo has been virialized. Otherwise, one should treat the 
recoil velocities in a more proper way as considered in 
Ref. ^3 t° estimate the resulting halo profile. Hence we 




FIG. 1: Solid (dashed) curves represent respectively from up 
to down the rotation curves (mass profiles) for the X halo 
which is from the NFW profile in the CCDM model, the Y 
halo in the DCDM model, and the case in which X decays 
into relativistic particles. The x-axis (y-axis) is in unit of 
r2oo of the X halo (V200 for solid curves and M200 for dashed 
curves) . 



can see that we have found out the allowed parameter 
space which is consistent with current observation data 
and justified the method we used. In short, the frac- 
tion of remaining X particles in the recently formed Y 
DM halo is tiny and given by £x ~ 10~ 9 , and they are 
decaying at t he p resent epoch to produce the observed 
UHECR flux [13. Furthermore, the possible distortions 
of the ionization history of the Universe caused by the en- 
ergy injection from decays of these relatively short-lived 
X particles have been recently discussed and the super- 
heavy DCDM model is able to provide a good fit to the 
current CMB anisotropy and polarization da ta ffil . On 
the other hand, the scenario proposed in Ref. [ll| , where 
£x = 1 and tx ~ to, would produce unacceptably large 
flux of UHECR unless the relativistic particles produced 
from the X decay involve some exotic quarks and leptons 
which are weakly interacting and may generate UHECR 
at an acceptable level by interacting with the interstellar 
medium when propagating to the Earth. 



IV. PHENOMENOLOGICAL IMPLICATIONS 

We have shown that the KR model that has attempted 
to explain the presence of UHECR with energies beyond 
the GZK cutoff can easily provide a DCDM solution for 
the problem of subgalactic structure formation in the 
CCDM model. In the DCDM model in which X DM 
decay into relativistic particles ^lj, not only halo core 
density is lowered but also small dwarf galaxies are dark- 
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ened due to core expansion and subsequent quenched star 
formation. It has also been argued that presently ob- 
served dwarf spheroidal galaxies with lower velocity dis- 
persions were resulted from decaying dark matter and 
subsequent core expansion in a small fraction of halos 
with high velocity dispersions [23|. This model predicts 
that the small-scale power at higher redshift is enhanced 
compared to the CDM model as well as the gas fraction 
in clusters should decrease with redshift. The latter can 
be tested by X-ray and Sunyaev-Zel'dovich effect obser- 
vations. However, this model has been criticized for that 
the reduction in the central density of clusters of galax- 
ies due to X DM evaporation might be too large to be 
compatible with observations and could even be harm- 
ful to the halo substructure formation 

El- 

It has been 

pointed out that this excessive reduction can be remedied 
if X particles decay into non-relativistic stable massive 
Y DM, and shown that the Y DM provides well fits to 
the rotation curves of low-mass galaxies and does not 
necessarily produce a significant reduction of the central 
DM density of certain dwarf spheroidals 0] ■ Undoubt- 
edly, detailed numerical simulations of the subgalactic 
structure formation in the DCDM model versus high- 
quality observations on the properties of subhalos and X- 
ray/Sunyaev-Zel'dovich effect of clusters would test the 
DCDM model and should differentiate the two scenar- 
ios. Remarkably, the subhalo astrophysics at kpc scales 
may provide a hint to understand the mass difference be- 
tween X and Y in the KR model at energy scale of grand 
unification. 

To test models of superheavy particles directly in ter- 
restrial particle accelerators is quite impossible. How- 
ever, the particle spectra and the arrival directions of 
UHECR produced from decays of superheavy particles 
in the Galactic halo can provide crucial tests. Super- 
heavy particles decay into ultra high-energy quark and 
lepton jets which fragment predominantly into photons 
with a small admixture of protons [lsl Ma. Although 
UHECR observations seem to show a subdominant pho- 
ton flux the photon flux with energies near the 
GZK cutoff may be attenuated in the cascading of the 
jets in the radio background and intergalactic magnetic 
fields [22| . The ultra high-energy neutrino flux accompa- 
nying the UHECR has been calculated [HE El to be 
much higher than the proton flux due to the long mean 
free path and high multiplicity of neutrinos produced in 
high-energy hadronic jets. This neutrino flux is near the 
detection limit of the on-going AMANDA neutrino exper- 
iment and will be severely constrained by the upgraded 
AMANDA and next generation neutrino telescope Ice- 
Cube. Because of the off-center location of the Solar 
system in the Galactic halo, some amount of anisotropy 
in the arrival directions of UHECR is expected [22| • Re- 
cently it was claimed that no significant deviation from 
isotropy is found, based on the data from the SUGAR 
and the AGASA experiments taken a 10-year period with 
nearly uniform sky coverage |33| . This may be overturned 
due to insufficient statistics. It is likely that the signal 



of the predicted anisotropy will have to wait to be tested 
by the upcoming Pierre Auger Observatory. 

As pointed out by KR [l!j, instanton mediated de- 
cay processes typically lead to multiparticle final states. 
Thus X particle decays will produce a relatively large 
number of quark jets with a fairly flat energy distribution 
and rather hard leptons as compared to typical perturba- 
tive decays of superheavy particles. This may leave a dis- 
tinct signature in the predicted UHECR spectrum which 
may help in distinguishing the KR model from other 
DCDM models. Furthermore, in the KR model which 
has 6 < 1, the energy of the relativistic Y particle is 
about mx /2 and the flux of Y particles in the Solar vicin- 
ity is approximately given by nxRhalo/rx ~ 10~ 5 nx , 
where Rhaio ~ lOOkpc is the size of the Galactic halo. 
This flux is about two orders of magnitude lower than 
the local flux of typical halo DM which is estimated as 
nxvx ~ KU 3 nx (where vx is about 300 km/s). If the 
Y particle interacts weakly with ordinary matter, it may 
scatter with the target nucleus with mass mjv in a cy- 
rogenic detector and deposit a huge amount of energy of 
order mjv(l — S)~ 2 in the detector. This deposit energy 
is much larger than that of a typical halo DM particle 
which is about rriNV x . This may give a new thought 
to the direct detection of halo DM. Unfortunately, since 
the local number density of X is nx ~ (GeV/mx)cm -3 
strongly suppressed by the mass of X and X is weakly 
interacting, the direct search for halo X particles or the 
indirect search for high-energy neutrinos from decaying 
X particles captured in the Sun or the Earth in current 
experiments are elusive |34j . However, it is worth noting 
that the fluxes of A-induced high-energy neutrinos from 
the Sun and the Earth are expected to be similar, though 
they are relatively low, to those considered in a different 
context of annihilation of strongly interacting superheavy 
DM which are predominantly tau neutrinos with a flat 
energy spectrum of events at about few TeV |35| and 
distinguishable from the energy spectrum of high-energy 
neutrinos induced by neutralino DM [3r| . 



V. CONCLUSIONS AND DISCUSSIONS 

In conclusion we have discussed the implication of 
the Kuzmin-Rubakov's decaying superheavy dark matter 
model for generating cosmic rays with energies beyond 
the Greisen-Zatsepin-Kuzmin cutoff to the subgalactic 
structure formation of the Universe. The model involv- 
ing a new SU (2)x gauge interaction and two left-handed 
SU(2)x fermionic doublets X and Y can easily accom- 
modate decaying dark matter scenarios for solving the 
cuspy halo problem inherent in the collisionlcss cold dark 
matter model. Intriguingly, the longevity of X particles 
due to instanton mediated decays may explain the pres- 
ence of a small cosmological constant as well. 

The drawback is that we require the near-degeneracy 
of X and Y particle masses. However, this may have 
recourse to physics at the relevant high energy scale. In 
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order to obtain the near mass degeneracy between X and 
Y particles, we assume that at high energies there is a 
symmetry, for example an exchange symmetry between 
X and Y, that makes their masses equal. Small mass 
differences could be generated by radiative corrections 
from symmetry breaking terms arising via threshold cor- 
rections near grand unification scale or even from stringy 
effects near Planck scale. For example, consider a term 
XiCi which contains X and other heavy fields. The one- 
loop correction lifts X mass by a factor of \\ /l€m 2 mx , 
giving rise to 8 ~ ICC 2 for Ax ~ 1. To get an even smaller 
6, we may use the idea of collective breaking of symme- 
tries. Instead of using one single coupling to break the 
symmetry, we introduce another similar coupling X2C2 in 
such a way that each coupling by itself preserves sufficient 
amount of symmetry such that the mass degeneracy be- 
tween X and Y is exact at one-loop level. It is only when 
the simultaneous presence of both symmetry breaking 
terms the mass degeneracy will be lifted. Therefore the 
radiative corrections which lift the mass degeneracy of 
X and Y are necessarily proportional to both Ai and A2. 
Hence this mass degeneracy splitting effect occurs at two- 
loop level and is of order Xf/Wir 2 X^/Wtt 2 which is suffi- 
ciently small even for Ai ~ A2 ~ 1. An alternative mech- 



anism for generating a small mass difference between X 
and Y particles is closely related to the result of instan- 
ton effects considered here. The mass relation between X 
and Y may be slightly modified by nonperturbative mass 
renormalization due to instanton- induced counterterms, 
similar to instanton-generated quark masses considered 
in QCD physics |37j . 

It is quite interesting to link different astrophysical and 
cosmological problems in a single particle model at grand 
unification scale. Future observations of dark matter ha- 
los and ultra high-energy cosmic rays, halo dark matter 
experimental search, and future CMB anisotropy and po- 
larization measurements will test the decaying dark mat- 
ter models and shed light on the mass degeneracy of X 
and Y. 
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